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We carry out a numerical study of the quantum Hall ferromagnetism in a two-subband system 
using a set of experimental parameters in a recently experiment [X. C. Zhang, I. Martin, and H. W. 
Jiang, Phys. Rev. B 74, 073301 (2006)]. Employing the self-consistence local density approximation 
, for growth direction wave function and the Hartree-Fock theory for the pseudospin anisotropy energy, 

' we are able to account for the easy-axis and easy-plane quantum Hall ferromagnetism observed at 

total filling factor v = i and v = 4:, respectively. Our study provides some insight of how the 
anisotropy energy, which highly depends upon the distribution of growth direction wave functions, 
determines the symmetry of the quantum Hall ferromagnets. 
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PACS numbers: 73.43.Nq, 71.30.-fh, 72.20.My 

I. INTRODUCTION 

Multi-component quantum Hall systems have exhibited a collection of interesting phenomena which are the man- 
ifestation of electronic correlations [1]. Such correlations become particularly prominent when two or more sets of 



I Landau levels (LLs) with different layer, subband, valley, spin, or Landau level indices are brought into degeneracy 
C [2, 3, 4, 5]. In experimental systems, different LLs can be tuned to cross by varying gate voltage, charge density, 
, ^ , magnetic field or the magnetic field tilted angle to the sample. One of the attractions is the formation of quantum 



Hall ferromagnets (QHFs) due to the exchange interactions of the two subbands states, termed as pseudospins [6, 7]. 
Self-consistent local density approximation (SCLDA) and Hartree-Fock mean field method can be performed on the 



> 

" calculation of ground state energy and quasi-particle energy gap [4, 7, 8, 9]. 

' when two LLs were brought into degeneracy. The QHFs can either be easy-axis or easy-plane, depending on the details 

OO : 

■ of LL crossing configurations [10, 11]. In this paper we follow the theoretical framework of Jungwirth and MacDonald 



Recent experiments in single quantum well with two-subband occupied systems [10, 11], showed evidence of QHFs 



. !^ ' [8] and numerically calculate the pseudospin anisotropy energy using the sample parameters in the experiment of 
5_j Zhang et al. [11]. The result confirms the QHFs taking place at total filling factor u = i and v ^ A are expected to 



be easy-plane and easy-axis QHFs, respectively. As we will discuss in this paper, the gate bias voltage, which affects 
the spatial distributions of the wave functions of the two subbands, plays a leading role in the formation of easy-plane 
and easy-axis QHFs. 
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II. PSEUDOSPIN QUANTUM HALL FERROMAGNETS 

The pseudospin representation is used to describe the valence LLs degenerated in the Fermi level [7, 8]. Following 
the theoretic studies of pseudospin QHFs [7, 8], here we focus on a two-subband two dimensional electron system, 
in which the LLs are labeled by s), where ^ = S/A is first/second subband (in the no biased quantum well, 
also called symmetry/antisymmetry subband) index, n — 0,1,--- is Landau level in-plane orbit radius quantum 
number, and s = ±i represents real spin. When two LLs are brought close to degeneracy but still sufficiently far 
from other LLs, one of them can be labeled as pseudospin up (cr ='{[) and the other as pseudospin down (cr =JJ.). In 
the experimental work of Zhang et al. [11], as shown in Fig. 1, we can label pseudospin up (cr =f|~) as (S*, 1, ■^) and 
pseudospin down (ct =JJ.) as (A, 0, i) at filling factor v — i. At filling factor = 4, we label pseudospin up [a ='f|-) as 
(S", 1, =F^) and pseudospin down (cr =-IJ.) as (A,0, ±^). Here the upper and lower signs before spin index refer to the 
crossing point at lower and higher magnetic field, respectively. 
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FIG. 1: Landau level diagram of two-subband system. Four Landau levels labeled by {£,, n, s) are crossing to each other during 
the increasing of the magnetic field (see text for details). Three arrows point out the crossing point of Landau levels at filling 
factor v — 3 and v — 4. 

At this moment, we would emphasize the essential similarity between two-subband system and ordinary bilayer sys- 
tem (double quantum well). On one hand, due to a separation of charges to opposite sides of the well originating from 
the Coulomb repulsion, a wide single quantum well can be modeled as an effective bilayer system [4, 8, 12]. Therefore 
in a single quantum well two-subband structure, the states of electrons can be characterized by two parameters similar 
to a bilayer system: the tunneling gap Ag^^s and effective layer separation d. The As^s is chosen as equal to the 
difference between the lowest two subband energy levels and d is given by the distance between two centers of the 
charge distribution in this single quantum well. On the other hand, let us assume the two dimensional electron gas is 
in the x — y plane in the following calculation and discussion, so the sample growth direction is aligned with the z axis. 
In a no biased bilayer system, although the layer wave function in z direction is spatially separated in two layers, once 
included in the effect of tunneling between two layers At and rediagonalized the z direction Hamiltonian, it will be just 




like the two-subband system which has symmetry and antisymmetry subbands. And the symmetry-antisymmetry gap 
^SAS is equal to the inter-layer tunneling Aj. Even when a bias voltage between two layers exists, the circumstance 
is more or less the same. Since they are theoretically identical, all the pseudospin language theory used in bilayer 
system is available in the two-subband system. 

While pseudospin up and pseudospin down LLs are degenerate but the number of electrons is not enough to fill all 
the two LLs, electrons will stay in a broken-symmetry ground state. Actually, the state electrons choosing is a linear 
combination of two pseudospin LLs which minimizes the system total energy. Typically, the many-body ground state 
can be written as follows [8]: 

i*N)=n4.jo), (1) 
fe=i 

where ct^ creates the single-particle state oriented in a certain unit vector m = (sin 6 cos ip, sin 6 sin ip, cos 6) with 
wave function: 
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i>rh,k{r) = COS {-^)tlJf,k{'^ + sin {-^)e"^tlJii.,k{r). (2) 

Here tpa,k{r) is the single-particle state wave function which contains growth direction subband wave function X^{z) 
and in-plane LL wave function Ln^s,k{x,y)- 

i}a,k{r) = \^{z)Ln,s,k{x,y). (3) 
Then the Hartree-Fock energy of the system can be obtained as following [8] : 

^..(rh) . = _ ^ U-\U.^.-\U.^^.^\ E (4) 

Whether the ground state is easy-axis or easy-plane only depends on the quadratic coefficient in pseudospin magne- 
tization rrii (i.e. the pseudospin anisotropy energy Uxx, Uyy and Uzz )'■ 

If (Pq 1 f dPq ,_. 

The first g = term in Eq. (5) is the Hartree term and the second term is the Fock term. Vij {q) can be expanded to 
sum of several pseudospin matrix elements WCTj.cr^.tri.cral'f)' which are products of subband and the in-plane parts [8]. 
If Uzz < Uxx — Uyy, thc systcm is in easy-axis QHF, which means the pseudospin magnetization m is aligned cither 
up or down; and if Uzz > Uxx = Uyy, the system is in easy-plane QHF, which is a coherent superposition of the two 
pseudospin LLs. 

III. GROWTH DIRECTION WAVE FUNCTION CALCULATION USING SCLDA 



We numerically calculate the growth direction (z direction) wave function A^(z) in Eq. (3) using SCLDA to compare 
the pseudospin anisotropy energy Uxx and Uzz in Eq. (4) [13]. In the SCLDA method, wave function X{z) is described 



by the Schrodinger equation: 
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+ Vfc(z) + Vgateiz) + V^ciz) + Vr (z) Ai(z) = £,Ai(z). 



(6) 



2m* dz^ 

Here m* is the effective electron mass in GaAs, VJ, corresponds to the conduction band discontinuity, Vgate is the 
bias potential caused by the difference of front and back gate voltage |AV^|, Vxc refers to the exchange-correlation 
potential related to the electron charge distribution n{z) (We use the form of Vxc given by Hedin and Lundqvist [14]). 
The Hartree term Vh due to electrostatic potential is given in the Poisson equation: 

2^e2 



Vniz) 



dz' \z — z'\ n(z'). 



(7) 



The subband energies e^, wave functions A^(2:) and the electron charge distributions n^i^z) of both subbands can 
be calculated by solving the Schrodinger equation Eq. (6) and the Poisson equation Eq. (7) simultaneously [15]. 



IV. NUMERICAL RESULT AND DISCUSSION 



Taking the unit of energy as /ds {Ib is the magnetic length), we calculate the pseudospin anisotropy energy 
Uxx, Uyy and Uzz at filling factor v = Z and v — A in Zhang et al.^s work [11]. At total filling factor v ~ 3, 
Uxx = Uyy ~ —0.370 < Uzz — 0.017, so the system will stay in easy-plane QHF. At one degenerate point of total filling 
factor 1^ = 4, Uxx = Uyy = > Uzz = —0.1266, and at the other degenerate point Uxx = Uyy = > Uzz = —0.1293. 
Both of the energy differences of Uzz — Uxx at = 4 indicate easy-axis QHF ground states. 
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FIG. 2; (Color online) Phase diagram of Uzz — Uxx (unit: e^/eZo)at filling factor (a) u = 3 and (b) u — A. In the Uzz — Uxx < 
(left/blue) region easy-axis QHF is energetically favorable while in the Uzz — Uxx > (right/red) region, easy-plane QHF is 
favorable. Black lines in each figure labeled by Uxx ~ Uzz show the critical positions where a quantum phase transition from 
easy-axis to easy-plane QHF occurs. The three stars in these two figures denote the v = ?> and = 4 experimental parameter 
in the experimental work of Zhang et al. 



To illustrate the evolution from easy-axis QHF to easy-plane QHF, we calculate the phase diagrams of Uzz — Uxx 
as a function of bias gate voltage [Al/^] and total density n at filling factor v = i (Fig. 2(a)) and v — A (Fig. 2(b)). 



In the following discussion, in order to make consistency, we choose e^/eZo ('o = lOnm) as the unit of energy. In 
the phase diagram (Fig. 2), we label the density and bias voltage |Ayf/| position where the crossing occurs in the 
experimental work of Zhang et al. [11]. In the left/blue (right/red) parts of each figure (Fig. 2(a) and Fig. 2(b)), 
where Uzz < {>)Uxx, easy-axis (easy-plane) has the lower electron interaction energy. From Fig. 2, we find that the 
anisotropy energy difference Uzz — Uxx is very sensitive to the bias voltage |Ay(/|. If we could vary the gate voltage 
across the black line labeled Uxx = Uzz from left to right in a determined density, a quantum phase transition from 
easy-axis to easy-plane QHF will happen. 
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FIG. 3; Energy components of anisotropy energy U^x and IJ^z v.s. bias voltage lAl/g]. Except for U^^ = at = 4, all the 
other terms are same at filling factor v = ?> and = 4. Uzzh and Uzzp correspond to the Hartree term and Fock term in the 
Hartree-Fock calculation of Uzz- 



In order to make a comparison of the effect of each term in anisotropy energy, we plot them in Eq. (5) as a function 
of |AF(7| at a certain density of 8.5 x 10^^ /cm? in Fig. 3. Here we have to point out that all the formulas at = 3 and 
4 are the same, except that Uxx vanishes a.t ly — 4. So the Fig. 3 is applicable to both filling factor v — i and v — A 
only if in the v — A situation the Uxx term is set to a constant zero. The reason for the different behavior of Uxx at 
filling factor v = i and v = A will be discussed in the following. Note that the dominated term in Uzz — Uxx is the 
Hartree term Uzzh, which is due to the electrostatic potential of electrons. It shows that the Uzzh term increases 
immediately with the increasing bias voltage jAV^gj. 

To examine what the role the bias voltage plays in the determination of easy-plane or easy-axis QHF, we give some 
SCLDA results in the Fig. 4. We define an effective subband separation Ad to describe the distance between the cores 
of first and second subbands wave functions. The \Ad\ as a function of bias voltage is plotted in Fig. 4, also with some 
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FIG. 4: (Color online) Effective subband separation lAdj as a function of bias gate voltage |AV^| at the density of 8.5 x 10^^ /cm^. 
Three insets show the well profiles (black dashed lines) with first (blue solid lines) and second (red dash-dot lines) subbands 
wave functions at different bias voltage: [Al/gl = O.Oy, 0.3^, and l.OV^. 



demonstrations of quantum well configurations and subbands wave functions in different bias voltages (|Al/g| — O.OV^, 

0. 3V^, and l.OF). The data in Fig. 4 are all selected from the same density 8.5 x 10^^/cm^ as in Fig. 3. It is obvious 
that the bias gate voltage changes the effective separation of two subbands while changing the well potential Vgate{z) 
in Eq. (6). The larger the bias voltage is added, the farther the lowest two subbands are separated. Since at a well 
separated z direction wave function configuration, all the electrons near the Fermi level filling in the same subband, 

1. e. the same pseudospin level, will raise a larger electrostatic energy, the easy-axis QHF is not favorable. Thus an 
easy-plane QHF can save more Hartree energy in the larger bias voltage situation. On the other hand, when the 
potential of the quantum well maintains a good symmetry in the small bias voltage limit, electrons can pick up one 
of the two pseudospin levels to keep the Hartree energy minimal and avoid the energetic penalty from inter-subband 
tunneling as well. 

As mentioned before, the difference of Uxx between total filling factor v — 3 and v — A could be explained in the 
similar way. The anisotropy energy Uxx or Uyy is constituted of Hartree part and Fock part (see Eq. (5)). In our 
numerical calculation, we find that the Hartree terms of Uxx and Uyy are always zero at both = 3 and v = A. And 
the only non-zero term in Uxx or Uyy is the Fock term Uxxp = UyyF < at filling factor v — i, which owes to the 
exchange interaction. It implies that at filling factor = 4, where pseudospin up a ='\\ (5,1,^^) and pseudospin 
down (T =JJ. {A^ 0, ±i) have opposite real spins, the easy-plane anisotropy, in which the electrons stay in both subband 
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equally, would cost much more exchange energy. But at v = 3 the pseudospin up a =f|- (5, 1, |) and pseudospin down 
a =J| {A,0, |) have the same spin. Then there is no such problem need to be considered. Therefore, the easy-axis 
QHF is more likely to happen at total filling factor v = 4 than v = 3. 

On the basis of above discussion, we may summarize as following. The bias gate voltage added to the sample 
changes the quantum well profile in the growth direction as well as the spacial separation of the lowest two subbands 
wave functions. For a larger bias gate voltage, the potential of the well is much skewer, so the two subbands wave 
functions locate in the opposite side of the quantum well (inset of Fig. 4). As a result, the Hartrce energy will arise if 
all the electrons stay in one narrow subband or pseudospin level. Thus a easy-plane QHF, in which the electrons fill 
the two pseudospin levels equally, is more energetically favorable at a large bias gate voltage. In addition, a state with 
opposite real spins will expend more exchange energy, so the easy-plane QHF is more easily to form at a pseudospin 
configuration in which the two pseudospin level have the same real spin (filling factor v = 3 in this paper). 

V. CONCLUSION 

Using self-consistent local density approximation method and Hartree-Fock mean field theory, we calculated wave 
function in the growth direction and the anisotropy energy in the two-subband quantum Hall system. The data shows 
great consistent with the observed easy-plane and easy-axis quantum hall ferromagnets at filling factor v = 3 and 
1/ = 4 in the experiment of Zhang et al. [11] . Also, by analyzing the numerical result, we give an easy-to-understand 
explanation about the anisotropy occuring at these filling factors. 
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